Reversal modes in magnetic nanotubes 
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The magnetic switching of ferromagnetic nanotubes is investigated as a function of their geometry. 
Two independent methods are used: Numerical simulations and analytical calculations. It is found 
that for long tubes the reversal of magnetization is achieved by two mechanism: The propagation 
of a transverse or a vortex domain wall depending on the internal and external radii of the tube. 



During the last decade, interesting properties of mag- 
netic nanowires have attracted great attention. Be- 
sides the interest in their basic properties, there is ev- 
idence that they can be used in the production of new 
devices. More recently magnetic nanotubes have been 
growrti^i^i* motivating a new research field. Magnetic 
measurements^ numerical simulations^ and analytical 
calculations 5 on such tubes have identified two main 
states: an in-plane magnetic ordering, namely the flux- 
closure vortex state, and a uniform state with all the 
magnetic moments pointing parallel to the axis of the 
tube. An important problem is to establish the way and 
conditions for reversing the orientation of the magne- 
tization. Although the reversal process is well known 
for ferromagnetic nanowires; 6 ' 7 ' 8 ' 9 ' 10 the equivalent phe- 
nomenon in nanotubes has been poorly explored so far 
in spite of some potential advantages over solid cylin- 
ders. Nanotubes exhibit a core-free magnetic configu- 
ration leading to uniform switching fields, guaranteeing 
reproducibility^ 5 - and due to their low density they can 
float in solutions making them suitable for applications 
in biotechnology (see [1] and refs. therein). 

Let us consider a ferromagnetic nanotube in a state 
with the magnetization M along the tube axis. A con- 
stant and uniform magnetic field is then imposed antipar- 
allel to M. After some delay time the magnetization 
reversal (MR) will start at any end. MR or magnetic 
switching can occur by means of different mechanisms, 
depending on the geometrical parameters of the tube. In 
this paper we will focus on the reversal process by means 
of two different but complementary approaches: numeri- 
cal simulations and analytical calculations. Their mutual 
agreement sustains the results reported in this study. 

Numerical Simulations. Geometrically, tubes are 
characterized by their external and internal radii, R and 
a respectively, and height, H. It is convenient to define 
the ratio (3 = a/R, so that (3 = represents a solid 
cylinder and (3 — > 1 correspond to a very narrow tube. 
The internal energy, E, of a nanotube with TV magnetic 
moments can be written as 

N N 
»=1 j>i 

where Eij is the dipolar energy given by Eij = 



with 



Tij the dis- 



tance between the magnetic moments \ii and , fi i the 
unit vector along the direction of /x^ and n,j the unit vec- 
tor along the direction that connects fi t and /x^ . Jy = J 
is the exchange coupling constant between nearest neigh- 
bors and Jij = otherwise. E a = — J2iLi Pi ' Ha is the 
contribution of the external magnetic field. In this paper 
we are interested in soft magnetic materials, in which case 
anisotropy can be safely neglected*^ In our simulations 
we consider Nickel tubes with \fi t \ = 0.61 /is, lattice pa- 
rameter ao = 3.52 A and J = 3.5 meV. Dimensions are 
R = 15 nm, length H = 0.5 /*m, assuming a growth along 
the [100] direction of a fee lattice. Tubes in the above 
mentioned range of sizes have at least 10 s atoms, and 
then numerical simulations at the atomic level are out of 
reach with present computational facilities. In order to 
reduce the number of interacting atoms, we make use of 
the scaling technique presented before ^ applied to the 
calculation of the phase diagram of cylindrical particles. 
Authors show that such diagram is equivalent to the one 
of a smaller particle with linear dimensions d! = dr^ be- 
ing x < 1 and 77 « 0.55, if the exchange constant has been 
also scaled as J' = \ J ■ It nas also been showni^ that the 
scaling relations can be used together with Monte Carlo 
(MC) simulations to obtain a general magnetic state of a 
nanoparticle. We use this idea starting from the desired 
value for the total number of interacting particles we can 
deal with, which based on the computational facilities 
currently available, is around 3000. With this in mind 
we have obtained x ~ 10~ 3 , that leads to nanotubes 
with around 10 3 atoms each. 

We will simulate the reversal process at temperature 
T = 300 K, using the scaling technique described above. 
MC simulations were carried out using Metropolis algo- 
rithm with local dynamics and single-spin flip methods. 13 
One interesting point to be considered is the effect of scal- 
ing on temperature at which the simulations are carried 
out. As explained by Bahiana et al^ to keep thermal ac- 
tivation process invariant under the scaling transforma- 
tion, the energy barriers must also be invariant, there- 
fore, temperature should scale as the volume, that is, 

T = x 3v T. 

We perform numerical simulations for tubes charac- 
terized by (3 ranging from 0.17 to 0.83, starting with the 
saturated magnetization along the cylindrical axis, the z 



2 



axis, with the external field H a applied in the — z direc- 
tion. Since the nucleation of a domain wall is more likely 
to occur at the ends, and we wish to follow the propa- 
gation of a single wall, 6 the field has not been applied 
to the last 12 nm of one end, following a pinning proce- 
dure used in experiments with microwires. For every (3 
we tried \H a \ = 1-3, 1.5 and 1.7 kOe, without any sig- 
nificant difference in the results. At least five seeds were 
used for each set of H a and f3. 

Continuum Model. Results of our simulations and 
previous results in wires^^^ show three main idealized 
types of MR, where M changes from one of its two en- 
ergy minima (M = M$i, with energy E F ) to the other 
(M = — Mqz, with energy E F ) by a path such that the 
energy barrier is the difference between the energy max- 
imum (-E max ) and the energy minimum. These mecha- 
nisms are illustrated in Fig. 1 and correspond to: Co- 
herent Rotation, C , where all the spins (local magnetic 
moments) rotate simultaneously; Vortex Wall, V, where 
spins rotate progressively via propagation of a vortex 
domain wall; and Transverse Wall, T, where spins ro- 
tate progressively via propagation of a transverse domain 
wall. 
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FIG. 1: Magnetic switching modes in nanotubes. 

To determine the magnetization reversal type for a 
given geometry, we calculate the energy for modes C , 
T and V finding the corresponding energy barriers. We 
adopt a simplified description of the system in which the 
discrete distribution of magnetic moments is replaced by 
a continuous one, defined by a function M(r) such that 
M(r)5u gives the total magnetic moment within the el- 
ement of volume 5v centered at r. The total energy for 
each mode is given by the sum of exchange and dipo- 
lar contributions, which now are taken from the well 
known continuum theory of ferromagnetism?^ The ex- 
change term is given by E ex = A J ^(Vmi) 2 dw, with 
rrii = Mi/Mo (i = x,y,z) the cartesian components of 
the magnetization normalized to the saturation value Mq, 
and A the stiffness constant. The dipolar contribution 
is Ed = (Mo/2) /M • VUdv, with U the magnetostatic 
potential^ We will proceed to describe the magnetiza- 
tion m fc and evaluate the total energy E k for each of the 
three modes defined above, with k — C, T or V. 
Coherent Rotation (C). For a uniformly magnetized 
tube with m z = cos 8, the exchange energy is zero, and 



the total energy corresponds to the dipolar one given by 
E c (9) = TTn M$HR 2 (l-l3 2 )[sm 2 6> + (3cos 2 0-l)N z ]f±, 
where the demagnetizing factor N z is given by^ 
2/H [°°\-e-i H 



N,. 



[J l {qR)^PJ l {l3qR)Ydq. 



J p represents a Bessel function of first kind and order p. 
Thus for long tubes^ E c (9) has minima at 9 = and 
9 = n. That is E c (0) = E c {n) = E F corresponds to 
the energy of a single domain ferromagnetic configuration 
along the axis of the tube. 

Vortex domain wall (V). In this mode the magneti- 
zation can be written as 

( z, < z < z w — w/2 

m v (z) = < m<t,(z)4> + m z (z)z, z w — w/2 < z < z w + w/2 
{ -i, z w + w/2 < z < H . 

In this expression z = z w denotes the position of the cen- 
ter of the domain wall - of size w - that is the height at 
which the magnetization lies totally in the xy plane de- 
scribing a perfect vortex. We model the magnetization 
inside the wall with the functional m z (z) — cos0(z), 
with 6(z) = 7r ((z — z w )/w + 1/2), which describes cor- 
rectly the instantaneous shape of the vortex wall in our 
numerical simulations. Using m x = —m^(z)sbx(j) and 
m y = m0 (z) cos <fi, the exchange energy results 



E 



it Aw ln(l//3) + tt 3 AR 2 (1 - p 2 )/w, 



which is independent of the position of the wall, z u 
dipolar contribution can be written as 

oo 

E v d = ir^M 2 R 2 J %[Jx{qR) - p.h{PqR)] 2 {f s 4 
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In the above equation, f s = f s (w, z w ) and /„ = f v (w, z w ) 
arc related to the surface and volumetric dipolar energies, 
respectively, and are given by 



fs 



fv 
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Transverse domain wall (T). Let us orientate the x 
axis along the transverse direction. Then, the magneti- 
zation inside the wall (see Fig. 1) can be written as, 

m T (z) = m x (z)'k + m z (z)z, z w — w/2 < z < z. w + w/2 , 

using the same geometrical notation as in the vortex wall, 
with w the width of the transverse wall. We model the 
axial component of the magnetization in this domain wall 
using m z (z) — cos9(z). The exchange energy results 



E 



ir 3 AR 2 (l- l3 2 )/w. 



The dipolar contribution gives 



with 



TTfl M ( 
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and g v {w,z w ) = f v (w,z w ). The dipolar volumetric c 
tribution of the V and T walls have the same functic 
form, although not necessarily the same value due to 
different wall widths of each mode, both denoted by 
With above expressions for the energies of the differ 
MR processes it is possible to obtain the reversal mo 
as a function of the geometry of the tubes. We si 
by calculating the energy maximum for each path, wf 
for the C mode is determined by 9 — 7r/2, whereas in 
V and T modes occurs when the wall is in the middle 
of the tube, that is for z w = H/2. The energy barrier 
AE k of each mode can be calculated as AE k = E^ ax — 
E F , with -E^ ax the energy of the A:- mode evaluated in 
their respective maximum. In cases V and T we have 
minimized the energy with regard to the wall width w, 
and then we search for the reversal mode which costs less 
energy. 

TABLE I: Results of simulations for different values of /3. 



p 


0.17 


0.33 


0.5 


0.67 


0.83 


Mode 


T 


T(V) 


V 


V 


V 



Table I shows the modes observed in our simulations 
for different j3 values. While a T mode is observed for 
[3 < 0.33, V appeared always for (3 > 0.5. For (3 = 0.33 
three seeds (of a total of fifteen) leads to a mix behav- 
ior: The reversal process start as T turning to V later 
on. The nucleation and propagation of the wall is mon- 
itored by the value of fhi(z) = Mi(z)/Mo, average value 
of the components of the magnetic moment at a height 
z, relative to saturation value. The position of the wall is 
determined by the maximum of (1 — \rh z \) and is shown 
in Fig. 2. In this figure we present snapshots at two 
different stages of the reversal process showing the prop- 
agation of the wall along the tube for (3 = 0.17 (a), 0.33 
(b) and 0.5 (c). The solid line represents the average ax- 
ial component of the magnetization (m z ) while the other 
two (in-plane) components are given by the dotted and 
dashed lines. When m x and m y both average to zero we 
face the V mechanism. When one or both of these com- 
ponents are non-zero, it is T that is observed. Case (a) 
shows a clear T behavior with a helicoidal rotation of the 
magnetization along the tube. Case (b) shows a T mode 
at stage ti, which is lost at ti where V is the only mode 
present. The last case, (c), illustrates a V mode along 
the tube. 
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FIG. 2: (Color online). Snapshots of the reversal on tubes 
defined by /3 = 0.17 (a), 0.33 (b) and 0.5 (c) at two different 
stages of the process. The abscissa represents the axial coor- 
dinate, z, along the tube and the ordinates gives the average 
components of the magnetization. fh x : dashed line (blue), 
rh y : dotted line (red), and m z : solid line (black). 
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FIG. 3: (Color online). Critical radius as a function of /3 for 
different magnetic materials. 



Now we turn our attention to the results provided by 
our theoretical model. Upon obtaining the lowest en- 
ergy barriers among the three modes we found that C is 
present only in very short tubes, namely, when H ps w 
or less. Since we are mainly interested in long tubes, 
as those used in experiments and applications mentioned 
above, we will focus our attention on T and V only. Be- 
cause it is the energy of the wall which determines the 
reversal mode, if H > w the height of the tube has no 
effect on the reversion process. In this case, our results 
show that the reversal mode depends on the internal and 
external radius, the stiffness constant A and the satu- 
rated magnetization Mq. For each (3 there exists a criti- 



4 



cal radius, R C (P), at which both energy barriers, T and 
V, are equal. For R < R c {(i) the tube reverses its magne- 
tization creating a transverse wall, while for R > R c ((3) a 
vortex wall appears. By equating both wall energies it is 
possible to obtain R c (j3), which is illustrated in Fig. 3 for 
Nickel (thick solid line), Permalloy (dashed line), Cobalt 
(dash-dotted line ) and Iron (thin solid line) . Parameters 
for Permalloy, Cobalt and Iron have been taken from ref. 
16, and for Nickel are the ones mentioned before. For 
R = 15 nm in Nickel we observe that for (3 < 0.33 a T 
mode occurs, a mixed reversal behavior is observed for 
/3 = 0.33 (denoted by a dot on the Ni line) and a V mode 
appears for (3 > 0.33. These results perfectly agree with 
our numerical simulations (see Fig. 2 and Table I). 




FIG. 4: (Color online). Solid lines (black) illustrate (a) the 
energy barriers and (b) wall width for the V mode and dashed 
lines (blue) correspond to the T mode. 



Fig. 4 illustrates our analytical results for the Ni tubes 
used in our simulations. Figure 4(a) depicts the energy 
barriers for the V and T modes as a function of j3. We ob- 
served that at /3 » 0.33 the reversal mode changes from 
T to V, in accordance to our previous numerical sim- 
ulations. Figure 4(b) illustrates the corresponding wall 
widths for modes T and V, which are of the order of the 
nanometers. Then, for very short tubes (H « w) it is 
not possible to accommodate a transverse or vortex wall, 
giving rise to a coherent mode of reversal. 

In conclusion, there are two basic reversion processes in 
long ferromagnetic nanotubes depending on the radii of 
the tube. The transverse mechanism is present in tubes 
with R < R C (P) while for tubes with R > R c ((3) it is 
always the vortex mode that controls the magnetization 
reversal. For values of R and j3 close to the boundary 
between the two phases, instabilities arise and changes 
between T and V modes can occur during the process of 
reversion. Because nanotubes with radius smaller than 
20 nm are difficult to fabricate at present,—^ we can 
conclude that the vortex reversal mode is the most ob- 
served. 
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